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LONG-TERM  GOALS 

Understand  morphology  changes  of  mud  tidal  fiats  and  provide  quantification  of  the  sensitivity  of 
these  changes  to  tidal  action,  river  discharge  and  shoreline  development. 

OBJECTIVES 

Demonstrate  the  use  of  a  community  numerical  model  for  estimating  the  tidal  fiat  morphology.  The 
ONR/Delft  Community  model  is  being  evaluated  as  a  physically-based  numerical  simulation  tool  for 
several  kinds  of  investigations.  The  objective  of  the  project  is  to  further  demonstrate  it  use  for  tidal 
flat/channel  systems. 

Investigate  the  relative  roles  of  tidal  action,  river  discharge  and  shoreline  development  on  morphology. 
From  a  model  tuning  perspective,  this  objective  includes  advancing  the  understanding  of  the  sensitivity 
of  the  model  to  parameter  value  adjustments  and  inclusion  or  exclusion  of  specific  sediment  transport 
processes  and  characterization  in  tidal  fiat  and  channel  systems. 

Integrate  model  analyses  with  observations  and  field  planning  efforts.  The  model  fidelity  will  be 
improved  by  incorporation  of  observational  data  in  terms  of  model  configuration,  the  assignment  of 
boundary  and  initial  conditions  and  for  calibration  and  validation  efforts.  There  is  also  a  benefit  to 
using  the  model  to  help  with  the  design  of  the  field  programs  and  the  interpretation  of  observational 
datasets.  This  has  been  demonstrated  in  other  studies  (e.g.  Hibler,  et  al,  2008). 

APPROACH 

A  three-phase  approach  will  be  used.  All  phases  are  model  configuration  and  evaluation  intensive.  In 
the  first  phase  the  model  was  configured  with  best  available  information.  Sensitivity  tests  were 
conducted  on  variations  on  discharge  to  circulation  in  the  Skagit  Bay.  In  the  second  phase. 


I 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1.  REPORT  DATE 

2QQg  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2008  to  00-00-2008 

4.  TITLE  AND  SUBTITLE 

Geomorphologic  Modeling  of  a  Macro-tidal  Embayment  with  Extensive 
Tidal  Elats:  Skagit  Bay,  WA 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Pacific  Northwest  National  Laboratory, 1529  West  Sequim  Bay 

Road, Sequim,WA, 98382 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  aS 

unclassified  unclassified  unclassified  Report  (SAR) 

5 

Standard  Form  298  (Rev.  8-98} 

Prescribed  by  ANSI  Std  Z39-18 


adjustments  to  the  model  will  be  made  to  intensify  foeus  in  portions  of  the  bay  where  field  aetivities 
have  oeeurred  and  where  they  are  planned  to  oeeur  in  the  future  (see  Figure  1).  Ineorporation  of  the 
observational  dataset  will  enable  further  model  improvement.  In  the  third  phase,  final  model 
evaluation  and  the  ineremental  phased  improvement  assessment  will  be  eompleted  and  reported  upon. 

This  work  will  build  upon  what  has  been  reported  for  numerieal  modeling  of  geomorphology  (Lesser, 
2004;  Tonnon,  2007;  van  Duin  et  al.  2004;  Van  Rijn  et  al  2007;  and  Mareiano  et  al.  2005).  The 
sensitivity  analysis  of  will  focus  on  bottom  type  classification  (sediment  type,  vegetation  coverage, 
etc.),  and  its  impact  on  sediment  transport  through  bottom  roughness  parameter  variation.  Process 
studies  on  idealized  channels  will  be  conducted  as  needed. 

Both  Mr.  Lyle  Hibler  and  Dr.  Adam  Maxwell  will  conduct  model  simulations.  Lyle  Hibler  will  focus 
on  project  management,  model  experimental  design,  circulation  modeling,  reporting  and  interaction 
with  other  program  participants.  Adam  Maxwell  will  focus  on  carrying  out  numerical  experiments, 
sediment  transport  and  morphology,  and  data  visualization  and  management. 

WORK  COMPLETED 

Work  completed  since  work  was  the  project  was  initiated  (May  1,  2008)  includes  development  of  the 
initial  mesh  (see  Figure  1)  and  incremental  refinements. guided  by  preliminary  results  and 
consideration  of  the  sampling  plans  provided  by  other  program  investigators  (Dr.  Jamie  MacMahan  at 
the  Naval  Postgraduate  School,  Dr.  Britt  Raubenheimer  at  Woods  Hole  Oceanographic  Institute,  and 
Dr.  Jim  Thomson  at  UW/APL).  Preliminary  simulations  focused  on  likely  conditions  during  August 
and  September  2008  for  coincident  field  programs. 

RESULTS 

We  have  begun  to  quantify  the  sensitivity  to  Skagit  River  discharge  rates  of  estimated  drifter 
trajectories  in  one  study  region  (in  Skagit  Bay  near  the  mouth  of  the  North  Fork  of  the  Skagit  River). 
We  have  learned  that  in  order  to  provide  the  most  relevant  estimates  to  the  program,  we  need  to  refine 
our  model  mesh  with  finer  resolution  near  the  field  studies  sites  (i.e.  to  begin  to  resolved  select  tidal 
channels).  We  anticipate  that  the  nominal  resolution  in  these  area  will  be  O(lOm).  In  order  to  achieve 
this,  we  will  need  to  use  the  model  options  including  mesh  decomposition  or  nesting  as  well  as 
incorporation  of  higher  resolution  bathymetry.  The  current  computational  mesh,  simulated  particle 
clouds  for  low  and  high  river  flow  rates  are  shown  in  Figure  1.  This  figure  also  includes  simulated 
drifter  tracks  at  possble  NPS  drifter  release  points  (planned  for  September,  2008),  UW/APL ’s  in  situ 
instrument  arrays  (planned  for  June,2009)  and  WHOI’s  CTD  and  ADV  (August,  2008)  field  sites. 
Figure  2  shows  the  simulated  drifter  tracks  and  particle  clouds  using  predicted  September  23,  2008 
tides  and  a  range  of  Skagit  River  flow  rates  (200  and  800  m^/s). 
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Figure  1.  Computational  mesh  developed  in  Phase  1  with  large  river  discharge  (pink)  and  low  river 
discharge  (black)  particle  clouds.  Locations  labeled  V  are  the  WHOI 8/2008 field  sites,  the  location 
labeled  SA  are  the  planned  UW/APL  6/2009 field  sites.  The  green  and  yellow  lines  are  simulated 
surface  drifter  tracks  place  in  region  close  (Craft  Island)  to  planned  drifter  release  by  NPS. 
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Figure  2.  Simulated  drifter  tracks  and  particle  clouds  for  high  (magenta)  and  low  (white)  Skagit 
River  flow  conditions.  WHOI  instrument  locations  are  shows  as  blue  triangles  and  UW/APL 

instrument  locations  shown  as  green  triangles. 

IMPACT/APPLICATIONS 

The  potential  future  impact  from  this  work  will  be  that  the  ONR/Delft  community  model  will  be 
further  evaluated  for  geomorphological  simulation  in  an  enviroment  that  is  of  interest  to  the  Office  of 
Naval  Reseach  and  in  DoD-Navy  where  the  software  is  already  being  used  for  other  applications. 

RELATED  PROJECTS 

Modeling  work  is  concurrently  being  done  to  support  a  program  entitled  Observations  and  Modeling 
for  Source  Characterization  (N000140810508)  with  Dr.  Mark  Moline 
at(  http;//www. marine. calpoly.edu/auv/)  and  in  collaboration  with  Dr.  Eric  Terrill  at  Scripps 
Institution  of  Oceanography,  and  Dr.  Ap  Van  Dongeren  at  WL|Deltares.  Lyle  Hibler  and  Adam 
Maxwell  will  be  using  the  same  numerical  modeling  software  to  investigate  the  circulations  Southern 
California  coast  with  focus  on  San  Diego  Bay  and  the  Tijuana  River  under  subcontract  to  California 
Polytechnic  State  University.  Integration  of  modeling  efforts  with  observational  datasets  will  be  the 
focus  of  this  effort. 
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